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ABSTRACT 

I derive a second-order local relation between the redshift- space mass density field and 
the real-space velocity field. This relation can be useful for comparisons between the 
cosmic density and peculiar velocity fields, for a number of reasons. First, relating the 
real-space velocity directly to the redshift-space density enables one to avoid the in- 
dependent reconstruction of the density field in real space. Secondly, the reconstruction 
of the three-dimensional velocity field in redshift space, questionable because of its 
vorticity, is also unnecessary. Finally, a similar relation between the galaxy density field 
and the velocity field offers a way to break the fi-bias degeneracy in density-velocity 
comparisons, when combined with an additional measurement of the redshift-space 
galaxy skewness. I derive the latter relation under the assumption of nonlinear but 
local bias; accounting for stochasticity of bias is left for further study. 
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1 INTRODUCTION 

Comparisons between the cosmic density and peculiar ve- 
locity fields of galaxies are potentially a powerful tool to 
measure the cosmological parameter Q. This is so because 
in the gravitational instability paradigm the density and the 
velocity fields are tightly related and the relation between 
them depends on Q. In the linear regime, i.e., when the 
density fluctuations are significantly smaller than unity, the 
density-velocity relation (DVR) is 



where b is the linear bias factor. Then the relation between 
the galaxy density and velocity field is 



<W(r) 



V • v(r) 



(3) 



S(r) = -f- 1 (Q)V-v(r). 



(1) 



Here, 8 is the mass density fluctuation field, v is the peculiar 
velocity field, /(fi) ~ Q, 6 and I express distances in units 
of kms -1 . Galaxies are 'test particles' which trace the ve- 
locity field induced by the mass distribution. However, there 
are both theoretical arguments (e.g., Kaiser 1984; Davis et 
al. 1985; Bardeen et al. 1986; Dekel & Silk 1986; Cen & Os- 
triker 1992; Kauffmann, Nusser & Steinmetz 1997; Blanton 
et al. 1998; Dekel & Lahav 1998) and observational evidence 
(e.g., Davis & Geller 1976; Dressier 1980; Giovanelli, Haynes 
& Chincarini 1986; Santiago & Strauss 1992; Loveday et 
al. 1996; Hermit et al. 1996; Guzzo et al. 1997; Giavalisco 
et al. 1998; Tegmark & Bromley 1998) that they are biased 
tracers of the mass distribution itself. In the simplest model 
of linear and local biasing, the galaxy density fluctuation 
field is 



Ogal 



bS, 



(2) 



Thus, a comparison between the fields in question based on 
linear theory yields not an estimate of SI directly but rather 
of a quantity (3 = /(f2)/6. This has been called the 'fi-bias 
degeneracy problem 'Fl 

To break this degeneracy in the IRAS— potent compar- 
ison, Dekel et al. (1993) used nonlinear effects. These effects 
in any case should be accounted for in analyzing the data, 
since the observed fields are currently smoothed over scales 
which are mildly nonlinear. Dekel et al. (1993) stressed that 
the limits on Q they derived 'are valid only in the approxi- 
mation that the linear biasing model is correct'; simply put, 
they assumed linear biasing. However, nonlinear corrections 
to the linear biasing relation (^) are of the same order as 
the corrections to the linear DVR (Q), so in general if one 
accounts for nonlinear dynamics one should also account for 



This terminology in the context of linear density— velocity com- 
parisons makes little sense to me, since in general we can speak 
about degeneracy of parameters only if we have a number of con- 
straints on them. This is the case when attempting to estimate Q 
from the multipoles of the redshift-space power spectrum. Here, 
we have just one constraint. Still, the terminology has become 
standard so I will use it. 
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nonlinear biasing as well. Such an analysis has been per- 
formed recently by Bernardeau et al. (1999; hereafter B99), 
resulting in a method of disentangling Q. from (a priori) non- 
linear bias. 

The method of B99 assumes that both fields are given in 
real space, while the observed density field, like, e.g., IRAS, 
is given in redshift space. In density-velocity comparisons, 
the real-space density field is commonly reconstructed from 
the redshift-space one using an iterative technique (Yahil et 
al. 1991, Strauss et al. 1992), though one-step procedures 
have been also invented (Fisher et al. 1995; Taylor & Valen- 
tine 1999). All these inversion methods aim at correcting the 
positions of galaxies for their peculiar velocities. However, 
to predict the velocities using an integral version of equa- 
tion (^) or its mildly nonlinear extension, one has to assume 
the value of (3. This is a serious drawback of real-space com- 
parisons: they necessarily assume the value of a parameter 
which is to be subsequently estimated. 

One might think that a way out of this problem is to 
perform the comparison in redshift space. Such an analysis 
has been done by Nusser & Davis (1994), but using linear 
theory only. Nusser & Davis (1994) applied the Zel'dovich 
approximation to model the redshift-space DVR, but subse- 
quently restricted their analysis to linear regime. The reason 
was, as they pointed out, that the redshift-space velocity 
field is rotational at second-order, so it cannot be recon- 
structed from its radial component only. Since in the present 
paper I aim at disentangling Q. and bias by using nonlinear 
effects in the DVR, the redshift-space analysis is also inap- 
propriate. 

A true way out of this dilemma is to relate the real-space 
velocity field directly to the redshift-space density field. Such 
a relation allows us to avoid at the same time problems with 
the reconstruction of the three-dimensional velocity field in 
redshift space and of the density field in real space. One 
might ask whether this makes any sense, since the fields in 
question are defined in different spaces. However, the point 
is that they are not: real-space and redshift-space are useful 
concepts, but in fact they are merely two different coordinate 
systems^. I will explain this in more detail later on. 

In the present paper I derive a second-order relation be- 
tween the redshift space density and the real space velocity. 
Next I demonstrate that, as in 'real-to-real' comparisons, 
second-order terms in this relation offer a method for dis- 
entangling Q, and bias. Now the method is self-consistent 
because it does not make any a priori assumption about the 
value of (3. I give completely worked out an example of non- 
linear but local bias. The case of stochastic bias is left for 
further study. 

The paper is organized as follows: in Section I derive 
the described DVR up to seco nd-o rder in perturbation the- 
ory. Specifically, in Subsection |2. l| I derive it in the case of 
unsmoothed fields; in Subsection |2.2| I extend it to the case 
of smoothed evolved fields, in a form directly applicable to 
so-called density-density comparisons. Next, I compute the 
expected scatter in the latter relation (Section ^). In Sec- 
tion ^ I derive a similar relation between the galaxy density 

t Obviously, the transformation from real space to redshift space 
does not conserve the volume element and in general may be even 
non-invertible. 



field and the velocity field under an assumption of a nonlin- 
ear but local bias. In Section ^ I show how such a relation can 
be used to break the f2-bias degeneracy in density-density 
comparisons. Summary and conclusions are in Section H. 



2 REDSHIFT-SPACE DENSITY VERSUS 
REAL-SPACE VELOCITY 

The purpose here is to derive a second-order local relation 
between the redshift-space density field and the real-space 
peculiar velocity field. Intrinsic motions of galaxies along 
the line of sight cause distortion in galaxy redshift surveys. 
On one hand, this distortion complicates measurements of 
the statistical properties of the large scale galaxy distribu- 
tion and also cosmic density-velocity comparisons. On the 
other hand, it contains valuable information on the structure 
of the peculiar velocity fields. Since the amplitude of pecu- 
liar velocities depends on £1, measurement of the degree of 
the distortion can serve as a method of measuring this pa- 
rameter (see Hamilton 1997: a review on linear redshift dis- 
tortions). Similarly to density-velocity comparisons, linear 
studies of redshift distortions yield merely an estimate of f3 
and in order to disentangle f2 from bias, nonlinear effects 
have been recently investigated (Taylor & Hamilton 1996; 
Heavens, Matarrese & Verde 1998; Scoccimarro et al. 1998). 

On large scales (<; several ft -1 Mpc), coherent flows 
lead to compression of structures along the line of sight. 
This compression reaches maximum at a turnaround ra- 
dius, where the structure appears totally collapsed along 
the line of sight and thereafter begins to invert itself. On 
small scales (^ a few h" 1 Mpc), high velocities of galaxies 
in clusters stretch out the structure in redshift space, cre- 
ating the so-called 'finger-of-God effects'. These effects are 
commonly interpreted as an evidence for virialisation of clus- 
ters. However, this is not necessarily the case. During three- 
dimensional highly-nonlinear gravitational collapse (but still 
before shell crossing), a structure attains infall velocities 
which are typically much bigger compared to its size multi- 
plied by the Hubble constant. Therefore, young unvirialized 
clusters, still partly in an infall phase, may also appear in 
redshift space as 'fingers-of-God' (Chodorowski 1990). 

A perturbative approach which I will adopt here is valid 
for scales bigger than turnaround. Then the 'compression 
effect' dominates and as we will see below, the redshift-space 
overdensity is increased compared to the corresponding real- 
space overdensity. 

2.1 Unsmoothed fields 

At the beginning my analysis will be similar to that of Heav- 
ens, Matarrese & Verde (1998), who studied nonlinear cor- 
rections to the redshift-space distortion of the power spec- 
trum. Let r and s be the coordinates in real and redshift 
space, respectively, with the Local Group at the origin. The 
redshift-space position of a galaxy is related to its real-space 
position by (Kaiser 1987) 



Here, u(r) = v r jr and velocities are expressed relative 
to the Local Group. Let p and p s be the density fields in real 
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and redshift space. The mean density may vary in space, due 
to selection effects; I will denote the mean densities as <f>(r) 
and 4>s(s). The density contrast in redshift space is then 
defined by 1 + <5 s (s) = p 3 (s)/(f> s (s); similarly for the real- 
space density contrast. 

Number conservation of galaxies implies p s (s)d 3 s = 
p(r)d 3 r. This yields 



[! + *.(«)] Ms) = 



u(r) 



[l + «'(r)] *x 



[l + 5(r)]<j>(r) : 



(5) 



where ' = 3/ dr. I will make the standard distant-observer 
approximation: that the survey is large enough so that the 
modes with wavenumbers k < r _1 are negligible. Then in 
the above expression we can neglect the term u(r)/r; fur- 
thermore, a Taylor expansion of <fi(s) yields <j>(r) plus a negli- 
gible correction. Second-order expansion of the above equa- 
tion is then 

5s(s) = 8(r) - u'{r) - S(r)u'(r) + [u'(r)] 2 . (6) 

I will now use the fact that in real space the density 
field is up to second order a local function of the two velocity 
scalars: the expansion (the velocity divergence), 8, and the 
shear, E. Specifically (Chodorowski 1997; cf. Gramann 1993, 
Catelan et al. 1995, Mancinelli & Yahil 1995) 



6(r)-- 
Here, 
E 2 = 

1j j 7 = 



' 2 E 2 (r). 



(7) 

(8) 
(9) 

e = v r -v, (10) 

and I have neglected the weak O-dependence. The symbol 
6ij denotes the Kronecker delta. Using expression ((?]) in ^ 
gives up to second order 



| (dvi/drj + dvj/dri) - \8l 



S s (s) 



-/ -1 0(r) - u'(r) + / -1 0(r)«'(r) + [u'(r)] 2 + 



■f- 2 [B 2 {r) 



E»] 



3rf 
2 



(11) 

Thus, the redshift-space density field at a position s is 
a function of the derivatives of the real-space velocity field 
at a corresponding position r. This is in fact a local relation 
between the redshift-space density and the real-space ve- 
locity. Real-space and redshift-space are useful notions, but 
actually they are merely two different mappings of the same 
physical space to 1Z 3 . The coordinate r and the correspond- 
ing coordinate s, given by equation (^), label the same point 
in physical space. Therefore, any local relation between the 
redshift-space quantity-one and the real-space quantity-two 
naturally expresses quantity-one at a position s in terms of 
quantity-two at a position r. However, this is inconvenient, 
because in the process of comparison, for each position r 
we should, given the velocity field, seek the corresponding 
s. Fortunately, the value of any real-space function g at r 
can be expressed as its value at f = s plus Taylor-expansion 
corrections, which will depend on the velocity. To second 
order we have 



g(r) = g{s) - u(s)g' '(«) . 



(12) 



This allows us to cast the RHS side of Equation ( |ll| ) , involv- 
ing the velocity evaluated at r, to another form, but also 



requiring only the velocity, evaluated at f = s. (I have used 
a 'check' mark to emphasize that r is a real-space point but 
in general different from r, which is related to s by eq. ^). 
The result is the relation between the redshift-space density 
and the real-space velocity, all evaluated at s: 

5 S = -r L o-u'+ [uir'e + u ')]'+£r a (e 2 - |s 2 ) .(13) 

What does it exactly mean? It means that if we select 
a point s in redshift space, to estimate the redshift-space 
density contrast at this point we need to evaluate the RHS of 
the above equation at a real-space point f = s. Conversely, 
if we select first a point r in real space, evaluating the RHS 
of the above equation at this point we estimate the redshift- 
space density contrast at a redshift-space point s = r. In 
other words, we can now treat the two fields as if they were 
given in the same coordinate system: in fact we are now 
comparing the density with the velocity at two (slightly) 
different physical points (since s/s), but the Taylor terms 
in the resulting DVR provide the necessary correction. 

I reiterate that equation ( |l3| ) involves real-space deriva- 
tives of the real-space velocity field. Having measured the 
peculiar velocity of a galaxy, we have a choice: we can as- 
sign the velocity either to the galaxy's distance, or to its 
redshift. The former procedure leads to construction of the 
velocity field in real space; the latter to construction of the 
velocity field in redshift space. Thus, the definition of the 
redshift-space velocity field, v 3 , is v 3 (s) = v(r). This field 
is not vorticity-free: V s x v 3 (s) — V s x v(r) 7^ V,- x v(r), 
the last quantity being the vorticity of the velocity field in 
real space, which vanishes. On the other hand, the difference 
between V s x v s (s) and V r x v(r), and hence the result- 
ing vorticity in redshift space, is already of second-order in 
velocity. In sum, while the linear velocity field in redshift 
space is irrotational the nonlinear one is not and therefore 
cannot be reconstructed from its radial component only. 

As stated earlier, in density-velocity comparisons a 
common procedure is to first reconstruct the real-space den- 
sity from the redshift-space one. Relation (|l^) provides a 
basis to avoid this step. The application of this relation 
to density-density comparisons is obvious: (a) take the ap- 
propriate derivatives of the real-space velocity field at r; 
(b) compare the result directly to the redshift-space den- 
sity contrast at s — r. As for velocity-velocity comparisons, 
relation |t3| ) can be straightforwardly used to reconstruct 
the real-space velocity field directly from the redshift-space 
density field. The real-space velocity field is irrotational thus 
describable as a gradient of some potential. Relation (|l^) re- 
duces then to a differential equation (in real space) for the 
velocity potential, with the source term S(r) = S s (s = r). 
This equation can be solved perturbatively: at linear order, 
by the method of Nusser & Davis (1994) ;0 at second order 
similarly, with the source term modified by the nonlinear 
terms in velocity, approximated by linear solutions. In the 
present paper I concentrate on the application of the rela- 
tion (|l^) to density-density comparisons; the application to 
velocity-velocity comparisons will be addressed elsewhere. 

The second-order DVR @ has been derived in a quite 
simple and natural way. There is, however, a more formal 



t The linear part of equation ( juj ) coincides with equation (6) of 
Nusser & Davis (1994) in the distant observer limit. 
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way of deriving it and, more importantly, generalizing for 
all orders. For convenience, introduce a scaled velocity field, 



(i + l/)r^ + ^(i + I/ + ^/ 2 )r 2 ^ 2 



(17) 



j—i 



(14) 



Equation (^|) then reads 

*.|*=r = [e + /«' + / [u{d + fa')] ' + ± (6 



')} 



(15) 



where for any function g, g = — / _1 <?- Applying the approach 
of Scoccimarro et al. (1998), in Appendix ^|l show that in 
general 



&s I s=r 



n 

n=0 



\r [u n (5+fu')] 



(n) 



(16) 



Here, the superscript '' n '' denotes l d n /dz n \ This equation 
will yield the desired DVR when we substitute for the real- 
space density contrast S(r) its local estimate from velocity; 
in other words, when we make use of the real-space DVR. 
Equation fttEj) is easily obtained from ( |l^ ) when the second- 
order relation (Q) is used and the series above is conse- 
quently truncated at second-order (n = 1) terms. However, 
extensions of real-space DVRs beyond second order have 
been studied both analytically and numerically (Nusser et 
al. 1991; Mancinelli et al. 1994; Chodorowski & Lokas 1997, 
hereafter CL97; Chodorowski et al. 1998; B99; Canon et 
al. 1999), resulting in a robust analytical description of the 
DVR in the mildly nonlinear regime.^] Specifically, B99 and 
Canon et al. (1999) found nonlinear formulas for density in 
terms of the velocity derivatives, which have proven to be 
very accurate, when tested against N-body simulations. One 
might thus think that using them in equation ([l6]) would 
yield a very good description of the mildly nonlinear local re- 
lation between the redshift-space density and the real-space 
velocity. 

In fact, this is not so simple and in the present pa- 
per I restrict the analysis to second order. The reason is 
that equations (|l(j) or ( |l5"| ) are not directly applicable to 
real data, since the observed density and velocity fields are 
always smoothed over some scale. As in real space, the ana- 
log of equation ( flB] ) - or of its higher-order extension - for 
smoothed fields will significantly differ from its unsmoothed 
counterpart (see the next subsection). Therefore, smooth- 
ing should be accounted for. This, however, turns out to 
be a non-trivial task, even in real space (CL97, Chodor- 
owski et al. 1998, B99), and calculations in redshift space 
are even more difficult. Restricting the analysis to second 
order will enable me to straightforwardly calculate the ef- 
fects of smoothing. Furthermore, as I will show it later, 
a second-order analysis is sufficient for the purpose of the 
present paper, i.e., outlying an idea of breaking the degen- 
eracy between Q and bias in density-velocity comparisons 
by nonlinear effects. Therefore, including orders higher than 
second would make the analysis more complicated than is 
necessary here. 

At the end of this subsection, let us compare contribu- 
tions of the linear and nonlinear terms in relation (|l3|). To 
simplify the comparison, consider a spherical top-hat over- 
density. Then E 2 = 0, 9 = const, it' = Id, and hence 



I.e., for density fluctuations of order a few. 



The real-space counterpart of the above 'redshift-to-real' re- 
lation can be obtained by setting the /(fi) factors in the 
parentheses equal to zero, yielding correctly a special case 
of equation ([j]): 



6 = -r 1 e+&r 



(18) 



The amplitude of the second-order term relative to the lin- 
ear one in this real-space relation is given for 9 = 1 by 
(4/21)/ -1 ~ 0.2/" 1 . Hence, while for a low-fi Universe 
real-space nonlinear effects can be significant, for a high-fi 
($7 ~ 1) Universe they are rather weak (Dekel et al. 1993). 
However, nonlinear effects in the relation (^) are signif- 
icantly stronger. Peculiar velocities increase the redshift- 
space overdensity at both orders, but while the linear term 
is increased only by 33% for Q = 1, the second-order term is 
increased by more than a factor of three. Thus, the relative 
importance of the nonlinear terms in the 'redshift-to-real' 
DVR is significantly greater than of those in the real-space 
DVR, and so is the sensitivity of the degeneracy-breaking 
method presented here. 



2.2 Smoothed fields 

As stated earlier, the density and the velocity fields recon- 
structed from observations are always smoothed over some 
(large) scale R. At linear order, the DVR for smoothed 
evolved fields follows immediately from its counterpart for 
unsmoothed fields. From the linear part of equation (|l^) and 
equation ( |lo| ) we have 5 R = V r ■ {i> R ) + f(u R )', because the 
operations of smoothing ( R ) and differentiation commute. 
However, the operations of smoothing and multiplication do 
not and the derivation of an analog of equation ( |l5| ) for 
smoothed fields including second-order terms is already a 
non-trivial task. Here I will derive the relation in a form 
applicable to so-called density-density comparisons, i.e., ex- 
pressing density in terms of velocity. In such comparisons, 
a large-scale velocity field is used to reconstruct the associ- 
ated mass density field, which is subsequently compared to 
the galaxian density field. To simplify the notation, I will 
from now on implicitly assume that all fields are smoothed 
at scale R and will drop the superscript R on all variables. 

Being a non-local operation, smoothing inevitably in- 
troduces some scatter in any local relation between non- 
linear variables. As an estimator of the smoothed density 
from the smoothed velocity we cannot therefore do better 
than to adopt the conditional mean, the mean redshift-space 
density contrast given in general a set of velocity derivatives, 
{^s)ls« 4 /ar -,S 2 o i /8r Sr ii , Second-order relation (Q for un- 
smoothed fields involves velocity derivatives up to second 
order, namely 9, E, u' (first order) and 9', u" (second). In 
principle, it would be then natural to compute the quantity 
(3s}\e £ a' e 1 a" ■ in practice, however, this is not the best 
thing to do. (Leaving aside extreme complexity of such a 
calculation!) The differentiation of the noisy peculiar ve- 
locity field should be done with great care and even so it 
still leads to dangerous biases (Dekel et al. 1998). The cur- 
rent incarnation of the potent algorithm, as an estimator 
of the real-space mass density from the velocity field uses 
a nonlinear formula but involving only first-order velocity 
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derivatives (Dekel et al. 1998, Sigad et al. 1998). Second- 
order derivatives of the velocity fields reconstructed from 
current catalogs of peculiar velocities probably cannot be 
reliably estimated at all. It would be thus more reasonable 

to compute the quantity (S s }\g j, ^/Jj 

However, a simple method of breaking the f2-bias 
degeneracy in density-density comparisons, which I will 
present later on, relies on a relation between 5 S and just 9. In 
such a relation, unlike in relation ( |l5| ) , both the (total) linear 
term and the second-order term will be a separable function 
of and the velocity derivatives.^ (See the simplified rela- 
tion [l?]). This will enable one to perform the comparison in 
a model-independent way, as described below. Therefore, as 
the only descriptor of the velocity field I finally choose the 
variable 9. An estimator of the redshift-space density con- 
structed exclusively from 9 has clearly bigger variance than 
an estimator constructed jointly from 9, S and u' . Still, this 
variance can be sufficiently reduced if the velocity field is 
sampled at a large-enough number of independent volumes. 

Taking the conditional mean of both sides of equa- 
tion @ we have 



(S s )\ § = f(u% + (9 + A 2 )\ § , 
where 

A a = /[«($ + /«')]' + £(fl a -§i 
Define now the variable 



7 



' + A 2 



(19) 



(20) 



(21) 



CL97 derived a general expression for the conditional mean 
for two mildly nonlinear smoothed variables, given Gaus- 
sian initial conditions, assumed also here. An additional as- 
sumption was that the variables are equal at linear order. 
Expanding 7 to second order we have 7 = 9\ + 9 2 + A2, 
where 9i are i-th order perturbative contributions to 9; in 
general is 0{9\). One can see that A2, like #2, is already 
0(9 2 ). Thus, the variables 7 and 9 are indeed equal at linear 
order. To second order for such variables (CL97), 



< 7 )| g = + a 2 



- (o 2 )) 



(22) 



Here, a 2 = (5*3 7 — S 3 §)/6 and {S 3 §) denotes the skewness 
of the variable 7 (9). Using the definition (^l|) of 7, at leading 
order we obtain: 



a 2 = ±(eiy 2 (elA 2 ). 

By isotropy we have 



(23) 



(24) 



This, together with equations (|2lJ) and (|22j), used in 
equation ( |l9| ) yields 

(Ss)\e = (1 + §/) $ + a 2 (9 2 - (9 2 )) , (25) 
or, transforming back to the plain velocity divergence, 
(S s )\o = - (1 + §/) r'O + a 2 f- 2 (9 2 - (9 2 )) . (26) 

" Chodorowski 1997 derived a related expression for the real- 
space density, (<5) \ g g ■ 

II Obviously, this is always the case when the velocity field is 
described by one variable. 



This is a second-order relation between the redshift-space 
mass density, at s = r, and the real-space velocity diver- 
gence, at r. The calculation of the coefficient a 2 for different 



smoothing windows is outlined in Appendix Bl. In analyzes 



of observational data the fields are commonly smoothed with 
a Gaussian window. The value of a 2 for a Gaussian smooth- 



ing function is given by equation (B13) 



The coefficient a 2 is a measure of the degree of non- 
linearity in the DVR. Since a method for breaking the fi- 
bias degeneracy, which I will present in Section 0, relies on 
the value of a 2 , we should compute it accurately. Smooth- 
ing clearly affects the value of a 2 (see Appendix |Bl[ ). For 
simplicity, con sider the case of a top-hat smoothing. From 
equation (B12) it follows that for = 1, the value of a 2 for 
unsmoothed fields, a 2 , overestimates the true value by al- 
most 40% for the spectral index n = — 1 and by more for 
bigger values of n. Thus, smoothing changes the value of the 
coefficient a 2 significantly so it was indeed necessary to ac- 
count for it. On the other hand, it is worth noting that the 
true value for SI = 1 and n = — 1 is still more than a factor of 



two bigger than its real-space counterpart, equation (B19) 



3 SCATTER IN THE RELATION 

A local relation between the density and the velocity di- 
vergence derived in previous Section has a scatter, present 
already at linear order. The r.m.s. value of the scatter is 
given by the square root of the conditional variance. At lin- 
ear order we have 



<(*.-<f.>li)% = f((u-\e) 2 ) 



y 2 ~9 2 -lf 2 9(u')\- e+ f 2 (u' 2 )\- e .{27) 



The quantity {u')\§ is given by equation (^|). The calcula- 
tion of the quantity {u 2 )\g is presented in Appendix^. The 
result is given by equation (JC5)) , where crj is the variance of 
the field 9. This yields ((S s - {S s }\ § ) 2 ) |- = -|/Vf, or 



((Ss - (6 s )\§f) 



1/2 



3V5 



f<?8 ■ 



(28) 



The scatter is of the order of the redshift-space correction 
to the linear term in equation (|25[) . When compared to the 
total linear term in equation (Eq), the scatter is relatively 
small: for 9 ~ erg, it is ~ 20% for Q = 1 and accordingly 
less for smaller £1. This implies that redshift-space density 
is well (though not as well as real-space density) correlated 
with real-space velocity-divergence. 



4 GALAXIAN DENSITY VERSUS VELOCITY 

As already stated in the introduction, it is currently clear 
that galaxies are biased tracers of the mass distribution. In 
this section I will derive a relation between the redshift- 
space galaxy density field and the real-space velocity field 
under an assumption of a nonlinear but local bias. This is 
only a toy model for bias because there are good reasons 
to believe that bias is in fact somewhat stochastic (Dekel & 
Lahav 1998; Pen 1998; Tegmark & Peebles 1998; Blanton 
et al. 1998; Tegmark & Bromley 1998). However, a number 
of important conclusions can be drawn already from this 
model. 
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Equation ( p^[ ) has been derived from conservation of 
galaxy numbers, so the densities appearing in it are in fact 
the galaxy densities <5s S ' and <5i s ' (redshift- and real-space 
galaxy density contrasts respectively) . Hence, a more general 
form of equation ([l6]) is 



oo 



(") 



(29) 



the form (0) implicitly assumes no bias between the distri- 
bution of galaxies and mass. In a local bias model, the real- 
space galaxy density contrast is assumed to be in general 
a nonlinear, but local function of the mass density contrast 
(Fry & Gaztanaga 1993; see also Juszkiewicz et al. 1995), 



<5 (s) (r) =Af[S{r)] 



(30) 



Given a (nonlinear) local formula for the real-space mass 
density in terms of the real-space (first-order) velocity 
derivatives, 8(r) = T [dvi/ drj{r)\, equation (129) yields 



OO 

]T ±f n {«" [M o Hdvi/drj) + fu'] } 



(n) 



.(31) 



r(g)| 



Up to second order, T is given by expression (^) and 

N [6{r)] = bS(r) + \b 2 (S 2 (r) - (S 2 )) . (32) 

Here, b and b 2 are respectively the linear and nonlinear 
(second-order) bias parameters. The term (<5 2 } ensures that 
the mean value of 5^ vanishes, as required. Using this 
in equation (|nj and truncating the series at second-order 
terms we obtain 

= {be + fu' + f [u{w + fu')]' 

+ ± b (P-lt 2 )+±b 2 (P-{P))}\ r . (33) 

This is a second-order local relation between the 
redshift-space galaxian density and the real-space veloc- 
ity, derived under the assumption of local bias. This re- 
lation is for unsmoothed fields. To account for the effects 
of smoothing and to enable a model-independent density- 
velocity comparison we proceed analogously to when con- 
structing the corresponding relation for the redshift-space 
mass density. Namely, we compute the expectation value of 
the redshift-space galaxian density given the velocity diver- 
gence. The result is 



(Si e) )\s = (&+§/) 



: + 4 S) 



(34) 



or, transforming back to the plain velocity divergence, 
(5^)\ e = -(b + y)f- 1 e + a^f- 2 (e 2 -(9 2 )) . (35) 



The calculation of the coefficient 



,(s) 



is outlined in Ap- 
pendix B2 . The result for Gaussian-smoothed fields is given 



by formula (B27). The coefficient depends on three pa- 



rameters: Q, b and b 2 . 

In section I argued that in the related mass-density 
versus velocity-divergence relation the nonlinear corrective 
term has fairly big amplitude. However, in the above expres- 
sion for t he galaxy density it may even vanish, if in equa- 
tion ( B27 ) for the coefficient a^ the nonlinear biasing term 
of negative sign cancels exactly the terms due to nonlinear 
dynamics. Negative nonlinear biasing of the IRAS density 



field has been inferred from a preliminary analysis of its bis- 
pectrum (Scoccimarro, private communication). This makes 
physical sense, since the IRAS galaxies are underrepresented 
in the cores of rich galaxy clusters, which is effectively a non- 
linear antibiasing operation on the mass density field. This 
may partly explain why in the IRAS— potent comparison 
the measured nonlinear effects are weak (Dekel et al. 1993), 
although in real-space analyses they are generally weaker 
anyway. 



5 BREAKING THE fl-BIAS DEGENERACY 

Second-order equation ( |3^ ) can be used to estimate fi sep- 
arately from bias in density-density comparisons. As first 
pointed out by Dekel et al. (1993), the key point of the idea 
of breaking the S7-bias degeneracy by nonlinear effects is not 
to correct for them in the density reconstruction, but rather 
to use them as additional information. The degeneracy- 
breaking method presented here is similar to that proposed 
by B99. However, the method is now self-consistent, because 
it no longer needs to make any a priori assumptions about 
the values of estimated parameters. To apply the method 
of B99, we should first reconstruct the nonlinear real-space 
density field. In order to do this, we would have to assume 
simultaneously the values of £1, linear bias and nonlinear 
(quadratic) bias. Here, we relate the real-space velocity field 
directly to the redshift-space density field and this problem 
disappears. 

First steps of the degeneracy-breaking analysis can be 
outlined as follows: the output of a POTENT-like reconstruc- 
tion machinery should be real-space plain velocity diver- 
gence (6), a quantity derived directly from the data. Next, 
the velocity divergence should be compared to the redshift- 
space galaxian density, also a model-independent quantity. 
Since Q and bias factors appear only in the coefficients of 
the density versus velocity-divergence relation (eq. |35|), this 
comparison is clearly model-independent. Specifically, the 
redshift-space galaxian density should be plotted against the 
velocity divergence and a second-order polynomial should 
be fitted to the (hopefully) apparent correlation. The fitted 
linear and quadratic coefficients give us two equations for 
three variables Q, b and b 2 . We need therefore an additional 
constraint on these variables. 

As this constraint we can adopt the large-scale galaxian 
density skewness, since it is a second order effect, involving b 
and b 2 . Galaxian skewness can be measured only in redshift 
space. Therefore, we need a theoretical relation between the 
redshift-space galaxy density skewness Sgf' (which we can 
measure) and the redshift-space mass density skewness S^s 
(which we can compute) . Such a relation will depend on the 
real-space bias factors b and b 2 . This dependence may be 
quite complicated, since the operations of bias and redshift- 
space mapping do not commute. 

The issue of non-commutivity of bias and redshift- 
space mapping was addressed in some detail by Scoccimarro, 
Couchman & Frieman (1998). Specifically, they compared 
the redshift-space hierarchical amplitude for the galaxy field, 
Qi s ^ , to that predicted if bias and the mapping commuted. 
(The hierarchical amplitude is defined as the ratio of the 
bispectrum to the relevant sum of products of the power 
spectra.) In the local bias model, the real-space galaxian hi- 
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erarchical amplitude, Q^ s \ is related to the real-space mass 
hierarchical amplitude, Q, by Q (g) = Q/b+b 2 /b 2 (Fry 1994). 
If the assumption of commutivity were correct, this would 
imply 



Q 



(g) 



,/6 + b 2 /b 2 



(36) 



where Q s is the redshift-space mass hierarchical amplitude. 



Scoccimarro et al. (1998) showed that the actual Q 



(g) 

WO) B11UWCU L11CIL L11C CLU L 110,1 

differs from Q s /b + b 2 /b 2 . However, the differences are quite 
small: at most (for some values of angle) 20% for Q. = 1 and 
respectively smaller for smaller Q. Moreover, as a function 
of angle, they fluctuate in sign. The galaxy skewness is es- 
sentially an integral over the galaxy hierarchical amplitude. 
Therefore, the differences between the actual skewness and 
its approximate value from an analog of equation (^) will 
largely cancel out. Using equation (j36|), a simple calculation 
yields 

S (g> = Ssi+g^. 



b 2 



(37) 



For the reasons stated above, we can expect this equation 
to be accurate to ~ 10% or less. 

It should be emphasized that approximate commutivity 
applies only to particular combinations of some statistical 
quantities and not to the quantities alone. In the case of 
skewness, this combination is the ratio of the third moment 
to the square of the second. It is straightforward to show that 
at leading order, the mass second moment is {8 2 } — (1+ §/+ 

\f 2 ){8 2 ) and the galaxy second moment is ([^s^] 2 ) = b 2 (l + 
|/? + ±/3 2 ){S 2 ), where /3 = f{Q)/b. Therefore, <$ g) ] 2 } / 
b 2 {8 2 ) and the difference for, say, 6 = 2 and Q, = 1 is « 35%. 

The redshift-space mass density skewness was computed 
for scale- free power spectra by Hivon et al. (1995). They 
found that it depends very weakly on Q and can be well 
approximated by the real-space skewness, except in the case 
that Q is close to unity and the spectral index n <; 0. For 
more realistic power spectra the calculation can be done in 
an analogous way and the result is also expected to be well 
approximated by the corresponding real-space value. Again, 
it means only that redshift-space mapping globally changes 
the third moment and the square of the second by nearly 
the same factor. 

The redshift-space galaxian density skewness of the 
IRAS density field was measured directly by Bouchet et 
al. (1993). Kim & Strauss (1998) pointed out that sparse 
sampling of the IRAS galaxies makes estimates of the skew- 
ness using the moments method biased low. However, they 
proposed a method to measure the skewness by fitting the 
Edgeworth expansion to the galaxy count probability distri- 
bution function around its maximum, properly accounting 
for the shot noise. Kim & Strauss showed using mock cata- 
logs that this estimate of skewness is robust to sparse sam- 
pling; they then used this method to measure the skewness 
of the IRAS density field. Hence, we know the value of the 
galaxian skewness if we are to perform the IRAS-potent 
comparison, or at least know an effective way to measure it. 

A simple form of equation (^) enables us to eliminate 
b 2 from equation (B27) for the coefficient a^' ^ (for Gaus- 
sian smoothing). Doing this and using the result in equa- 
tion (B5|) yields 



c 2 + 



0: 



-(g) 



/r 2 (e 2 



o-g) 



(38) 



where C2 is given by equation (B28). The parameters £1 and b 



' + 



appearing in the above relation are not degenerate: while in 
the linear coefficient they appear only via their combination 
(5, in the quadratic coefficient they clearly do not. Therefore, 
measuring the linear coefficient and the parabolic correction 
in the galaxian density versus velocity-divergence relation 
one can, at least in principle, measure Q and b separately. 



6 SUMMARY AND CONCLUSIONS 

The main purpose of the present paper is to outline a self- 
consistent method of breaking the fi-bias degeneracy in 
density-velocity comparisons by nonlinear effects. A prob- 
lem which was present in previous similar attempts (Dekel et 
al. 1993; see also CL97) is that they assumed linear biasing. 
However, if we account for nonlinear dynamics, in general 
we should account for nonlinear biasing as well. There is also 
another problem which plagues all real-space comparisons, 
even these performed at the linear level, aiming merely at 
inferring the estimate of /3. Namely, in the process of recon- 
struction of the real-space density from the redshift-space 
density one has to assume the value of /3, the parameter 
which is to be subsequently estimated. Fj 

With nonlinear biasing accounted for in the compari- 
son, the resulting system of equations for bias parameters 
and Q, is underconstrained. As recently pointed out by B99, 
an additional constraint on f2 and bias can be inferred from 
the density field alone (namely, its skewness), so additional 
types of observations are unnecessary. In other words, there 
is enough information in the density field and the corre- 
sponding velocity field to disentangle SI from bias. As for the 
problem with the real-space density reconstruction, Nusser 
& Davis (1994) proposed to solve it by performing the com- 
parison directly in redshift space. 

One might thus think that to break the fi-bias degen- 
eracy in density-velocity comparisons by nonlinear effects 
one should: (a) account for nonlinear biasing, (b) add an 
additional constraint on it from the density field (e.g., the 
skewness) and (c) perform the comparison in redshift space. 
However, here we face a new problem: the redshift-space ve- 
locity field is rotational at second order, so in the nonlinear 
regime it cannot be reconstructed from its radial compo- 
nent only. Since the method of breaking the degeneracy is 
based on nonlinear terms in the density-velocity relation, 
the redshift-space analysis is also inappropriate. 

A way out of this dilemma is to relate the redshift-space 
density field directly to the real-space velocity field. Such a 
relation allows one to avoid simultaneously problems with 
the reconstruction of the density in real space and of the 
three-dimensional velocity in redshift space. In the present 
paper I derived this relation explicitly up to second-order 
terms. The relation is local, i.e. it expresses density at a given 
point in terms of the velocity derivatives at the same point. 
(Strictly speaking, it expresses the redshift-space density at 

** When attempting to reconstruct the nonlinear real-space den- 
sity field, this problem would get even worse: one would have to 
assume simultaneously the values of /3, linear bias and nonlinear 
(quadratic) bias. 
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a point s = r in terms of the real-space velocity derivatives 
at r.) Therefore, it can be straightforwardly applied to so- 
called density-density comparisons and I did this here. It 
can be also applied to velocity-velocity comparisons; I will 
address this point elsewhere. 

Next, I rigorously accounted for the effects of smoothing 
of the evolved density and velocity fields. This was necessary 
because smoothing changes the degree of nonlinearity in the 
density-velocity relation considerably. 

Finally, I derived a second-order relation between the 
redshift-space galaxian density and the real-space veloc- 
ity divergence, under an assumption of nonlinear but local 
bias. This relation (eq. ^) enables one to perform density- 
velocity comparisons in a model-independent way. When 
combined with an additional measurement of the skewness 
of the galaxy density field (eq. it also enables in prin- 
ciple to solve for Q and bias separately in density-density 
comparisons. In short, the idea is to measure the amplitude 
of the nonlinear corrective term in the relation. 

The density versus the velocity-divergence relation has 
a scatter even at linear order. This scatter is relatively small 
when compared to the linear term in equation (B8|). This 
means that the linear estimate of j3 is affected by the scatter 
in a rather weak way. In any case, a fairly large number 
of independent volumes in the comparison will enable to 
suppress this scatter, allowing one to estimate the linear 
coefficient and hence j3 accurately. 

The idea of breaking the f2-bias degeneracy presented 
here is somewhat related to that given by Pen (1998). In 
order to measure bias and fi separately, Pen proposed to 
use simultaneously redshift distortions of the power spec- 
trum and of the skewness. However, while in his method 
the estimate of il is based mostly on the redshift distor- 
tion of the galaxy skewness (which is quite a small effect), 
here it is based on the nonlinear effects in the redshift-space 
density versus real-space velocity relation. These effects are 
considerably stronger than in the corresponding real-space 
relation. 

Quality of the current velocity data is probably too 
poor to perform a degeneracy-breaking analysis similar to 
proposed here. The scatter in the (real-space) density- 
velocity diagram of the IRAS-potent comparison by Sigad 
et al. (1988) is huge and it results in significant errors in the 
estimate of /3. If even the linear coefficient is to some extent 
uncertain we cannot hope to reliably estimate the nonlinear 
corrections. 

The analysis in its present form is still very simplistic. 
Points that should be addressed in further studies in the 
first place are: (a) stochastic nature of biasing; (b) higher 
than second-order effects, which may play a role for scales 
smaller than ~ 12 ft -1 Mpc; (c) relaxing the distant-observer 
assumption, not strictly valid for current density-velocity 
comparisons.^] Recently, Taruya & Soda (1998) have de- 
rived a mildly nonlinear galaxy density-mass density rela- 
tion in a stochastic bias model. The calculation of the galaxy 
density-velocity divergence relation in such a model is math- 



tt Nevertheless, Kaiser's classical formula for the redshift-space 
linear power spectrum, derived in the distant-observer approxi- 
mation, is commonly used in analyses of redshift distortions in 
galaxy surveys like IRAS. 



ematically entirely analogous. Higher-orders effects are com- 
plex, but tractable; here I computed them for unsmoothed 
fields. Similarly, calculation of the density-velocity relation 
without the simplifying assumption of a distant observer, 
though more complex, is certainly possible. 

The purpose of this paper was to perform a first step 
towards self-consistent breaking the fi-bias degeneracy in 
density-velocity comparisons by nonlinear effects. I showed 
that the effects in the comparisons proposed here are fairly 
strong and that the fact that they are of the same order as 
possible nonlinearities in the biasing relation is not a funda- 
mental obstacle. The method of breaking the f2-bias degen- 
eracy presented here can be worked out in full detail when 
data of sufficient quality exists. I mean primarily catalogs of 
less noisy, more densely and homogeneously sampled pecu- 
liar velocities. 
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APPENDIX A: DERIVATION OF EQUATION @ 



In this appendix I present the derivation of exact relation (|16[), expressing the redshift-space density at s = r in terms of 
the real-space density and velocity at r. In the distant-observer approximation, the line-of-sight is taken as a fixed direction, 
which will be denoted here by z = t%. The real-to-redshift mapping (Q) reduces then to s = r — fu(r)z, where u is the 
z-component of the scaled velocity field, defined by equation From the conservation of the number of galaxies we have 

1 + 5 s (s) = J(r)[l + S(r)], where J(r) is the Jacobian of the mapping. In the plane-parallel approximation, J(r) — 1 — fu'(r) 
(exactly), where ' = d/dz. This yields 

S(r) + 1- J(r) 5(r) + fu'(r) 

Ss(s) = — W) — = i-mr) (A1) 

(eq. 3 of Scoccimarro et al. 1998). Fourier transforming the above equation we obtain 

*.(fc) = J ^Le-* k °6 s (s) = J ^J_ e -*-r e m,*(r) [ 5(r) + K(r) ] (A2) 
(eq. 4 of Scoccimarro et al. 1998). Thus we have expressed the Fourier transform of the density contrast in redshift spa ce as 



a real-space integral involving the real-space density and velocity fields. The inverse Fourier transform of equation (A2) is 



8 3 (s) = J d 3 ke tks 6 s (k) = J ^0^- [8(r) + /fi'(r)] e lfk ^ (r) e lk - (s - r> , (A3) 
where I have changed the order of integration. Expanding the first exponent of the integrand we have 

e iftu(f) e »(«-r) _ [fu{r)] n v« ifc.(«-r) _ (-1)" f n u™ (r) d" ^.(s-r) 

nl { z> ^ n\ dz n ' { ' 

n=0 n=0 

hence 



d k ik-(s-r) 



e ^ a - r > . (A5) 



The integral over k yields the Dirac delta distribution, So(s — r). Integrating the remaining integral over r by parts we obtain 

oo 

5 ^ = Hh r ^ n{s + {n) L=. ' (A6) 

where the superscript denotes 'd n /dz n \ Thus, if we select a point s in redshift space, to estimate 5 s (s) we evaluate the 
expression on the RHS of the above equation at a real-space point f — s. Conversely, if we first fix a point r in real space, 
evaluating the RHS of the above equation at r we estimate the redshift-space density contrast at a redshift-space point s = r. 
This yields equation (|li|). 



APPENDIX B: CALCULATION OF THE COEFFICIENT a 2 
Bl The case of no bias 

Our purpose here is to compute the coefficient a2, given by equation (p3|), for smoothed evolved fields. The calculation is 
much easier when performed in Fourier space. Denote the Fourier transform of the linear real-space density field <5i by e(fe). 
According to linear theory Q\ = 5i, so e(k) is also the Fourier transform of the linear scaled velocity divergence. We have 

r)2\ 2/ 



{el) = * 2 (R) , 



(Bl) 



where cr 2 (R) is the linear variance of the real-space density field smoothed over a scale R. To compute the quantity {§fA2^, 
we need the Fourier transform of A2 (defined by eq. Eoh. Following Heavens et al. (1998) I find that 



A 2 (fc) = 
where 



d 3 pd 3 g 



8 D (k-p- q)T 2 (p, q)e(p)e(q) , 



T 2 (p,q) = ~ (l-cos 2 0) + L 



2 2 



l-lpfJ>q 



£ + 1 

q p 



+f 



2 2 



2P 



2<7 
P 



(B2) 



(B3) 



In the above, cos$ = p ■ q and for any vector k = k/k. The remaining quantities are fi p — p ■ s and fi q — q ■ s, where s 
is a unit vector from the observer to the galaxy; in the distant-observer approximation one assumes s is constant within the 
smoothing radius. 

For Gaussian initial fluctuations, straightforward calculation yields 
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i (6l A 2 ) = J ^0JLw(p)W{q)W{\p + q\)P(p)P(q)T 2 (p, q) (B4) 

(see e.g. Juszkiewicz, Bouchet & Colombi 1993). Here, P is the real-space linear mass power spectrum and W is the Fourier 
transform of the window function. The calculation of the above integral with the kernel T2 given by equation ( |B3[ ) is similar 
to the calculation of the skewness of the redshift-space density field. Following Hivon et al. (1995) I note that the angular 
integration of the above expression for a fixed # is, up to a multiplicative factor, equivalent to averaging over different 
orientations of s. We then have p — (0,0, l)p, q = (sin$, 0, cos -£))q and s = (cos (j>' sin sin <j>' sini?', cos #')s, so 

Hp = cos , (B5) 

fj, q = cos <j>' sin i9' sin $ + cos i9' cos , (B6) 

and the averaging over the polynomials of /x p and [i q given $ is effectively over the trigonometric functions of cj>' and . The 
result is 



(T 2 ( P , q))^ q U = l + L + £ + f L + £ ) cos J E + J ) - ( I- ^ , cos' , . (B7) 



Using equations (Bl) and (B4) in equation (|23j) yields 
fl2= / ^$ W/(p)M/(9)W ' (b + ,7l)P(p)P(9)<T2(P ' 9)) ' (B8) 



where (T2) is given by expression (B7). The calculation of the coefficient 02 is now entirely analogous to the calculation of 
the skewness of the real-space density field. This is so because the averaged kernel has a similar functional dependence on 
p and q to that of the second-order solution for the real-space density field (see e.g. Juszkiewicz et al. 1993), with different 
coefficients. In the case of unsmoothed fields W(k) = 1 and we obtain 

4 Un) = | + { + f (B9) 



Comparing (B9) with expressions jl^ ) and (pi|), we see that a^ n ^ is equal to the value predicted by the spherical collapse 
model. 

To study the effects of smoothing I will consider power spectra with a power-law form 

P(k)=Ck n , -3<n<l, (BIO) 

where C is a normalization constant. For more realistic spectra, the value of 02 can be very well approximated by the result 
for scale-free spectra with the effective index at a smoothing scale defined as 

R da 2 (R) ^ . 

naB = ~^^R-~ 3 (BU) 

(CL97). For the case of a spherical top-hat smoothing function, using the results of Bernardeau (1994) I find that 

,(TH) _ fun) , „x ( f f 



C» + 3)(£ + fe). (B12) 



For the most interesting case of Gaussian smoothing, using the results of Lokas et al. (1995) I obtain 

(B13) 



■'■ _!_,.-„.„./ l> (3/2 ) - |F(5/2) 1 ' /2 : '"'- S 



... - —F<5/2) r - 
2 21 w ; 6 



+ — 
10 



f (3/2) ^—F(5/2) 



Here, 

FW»-*(=±2.2±i.|.i). (BH) 
^fl,.^ (• + »,- + »,■.') ,B I5 , 

and 2F1 is the hypergeometric function. The hypergeometric function can be expanded in a series of powers of n + 3 which 
converges rapidly. Fitting the expansion by a polynomial in n + 3 in the range —3 < n < 1 yields low-order accurate fits. 
Namely, 

F(r) = 1 + 4 r) (n + 3) 2 + d { {\n + 3) 3 , (B16) 
where 

4 3/2) = 0.03239, df 72 ' = 0.009364 (B17) 
and 

4 5/2) = 0.02183, 4 5/2) = 0.003531 . (B18) 
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The fit to the function F(3/2) is accurate to 0.6% or less in the range —3 < n < 1; the fit to F(5/2) is accurate to 0.2% or 
less. 

As in the case of top-hat smoothing, for n = —3 we have — a^" 1 ' and for other values of the spectral index 

. Thus, only for n = —3 are the values of the coefficient a 2 for smoothed fields equal to that given by the spherical 
collapse model. 

Finally, we recover the real-space values of a% (CL97, B99) by setting / = in equations (|B9|), (B12) and (B13): 



(TH) 



2, real 2^ 

and 



(I 



(G) 



-F(5/2) . 



(B19) 
(B20) 



B2 The case of bias 



Our purpose here is to compute the coefficient , entering in the relation ( |34j ) between the redshift-space galaxian density 
and the real-space velocity divergence. All calculations are similar to those in the case of no bias, presented in the previous 
Subsection. Using equation (63), analogously to equations (JlM-pCT) we can write 



(*(<% = 6 (l+f/3)0 + <A^)| 
where 



As) 



f [u(w + fa')] ' + ±b (e 2 - |e 2 ) + \b 2 (e 2 (e 2 )) 



The coefficient af is given by the following equation (cf. eq. |23j) : 



(B21) 
(B22) 



,(g) 



2 '0?AW 



The Fourier transform of a! 8 ^ is 



(B23) 



' [ ' ! ' ! ' ]iq 5 D (k-p-q)T 2 (p,q)e(p)e(q) - 4ty 3 b 2 a 2 8 D (k) . 



(2tt) 2 

where already an averaged form of the kernel is 

26 . 36 2 bf f 2 ( bf f 2 



7 



+ 



15 



6 + 10 



i P . Q 

COS 17 I 1 



2I> 
T 



2f_ 
15 



In the expression ( B24), th e last term is the Fourier transform of the term (fe/2) (6 ). 
Using equations (B2J) and (B24) yields (cf. eq. B8) 



4 S) = / ^0W(p)W( q )W(\p + q \)P(p)P( q )(T^(p, q ))- b -^. 
For Gaussian-smoothed fields, the result of the integration is 



(g,G) b 2 ,2 

a 2 2 2 ' 

where 



F(3/2) - -F(5/2) 



+ 



10 



F(3/2) 



3n + 8 



F(5/2) 



COS 1? . 



(B24) 



(B25) 



(B26) 



(B27) 



(B28) 



= f(Q,)/b and F(3/2) and F(5/2) are defined in equations (B14j)-(B15). In the limit of no bias (6 = 1 and 6 2 = 0), 



equations (B27)-(p28j) reduce to equation (D313J), as expected. The result for a top-hat smoothing can be obtained by replacing 



the factors F(3/2) and F(5/2) by unity. The result for unsmoothed fields then follows from setting the spectral index equal 
to -3. 



APPENDIX C: GAUSSIAN CONDITIONAL AVERAGES 

Our purpose in this Section is to compute the quantity (u' 2 }\§ at leading, i.e. linear, order. We may assume therefore that v! 
and 6 are Gaussian variables. For such variables, 

{u% = r^9. (CI) 
Here, a\, (ct|) is the variance of the variable v! (6) and r is the correlation coefficient defined as 
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r = 



Towards breaking the Omega-bias degeneracy in density-velocity comparisons 

(u'0) 



Using results of Appendix Bl we have a\i = (1/5)<t| and r = \/5/3. This yields 



a result which is obvious by isotropy. The conditional variance is 
(u' 2 )\ g - (w')|| = (l-r 2 )cr?, 
hence finally 



J2\ J2. _ _£ 2 

45 CT e ■ 
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